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Abstract 

Recent progress in understanding the physics of magnetic reconnection is conveniently summa- 
rized in terms of a phase diagram which organizes the essential dynamics for a wide variety of 
applications in heliophysics, laboratory and astrophysics. The two key dimcnsionless parameters 
are the Lundquist number and the macrosopic system size in units of the ion sound gyroradius. 
In addition to the conventional single X-line coUisional and collisionless phases, multiple X-line 
reconnection phases arise due to the presence of the plasmoid instability either in coUisional and 
collisionless current sheets. In particular, there exists a unique phase termed "multiple X-line 
hybrid phase" where a hierarchy of coUisional islands or plasmoids is terminated by a collision- 
less current sheet, resulting in a rapid coupling between the macroscopic and kinetic scales and a 
mixture of coUisional and collisionless dynamics. The new phases involving multiple X-lines and 
collisionless physics may be important for the emerging applications of magnetic reconnection to 
accelerate charged particles beyond their thermal speeds. A large number of heliophysical and 
astrophysical plasmas are surveyed and grouped in the phase diagram: Earth's magnetosphere, 
solar plasmas (chromosphere, corona, wind and tachocline), galactic plasmas (molecular clouds, 
interstellar media, accretion disks and their coronae. Crab nebula, Sgr A*, gamma ray bursts, 
magnetars), extragalactic plasmas (Active Galactic Nuclei disks and their coronae, galaxy clus- 
ters, radio lobes, and extragalactic jets). Significance of laboratory experiments, including a next 
generation reconnection experiment, is also discussed. 

PACS numbers: 



I. COLLISIONAL AND COLLISIONLESS RECONNECTION 



It has been a long held view that magnetic reconnection is primarily characterized by 
plasma coUisionality. This is evidenced by the common uses of the resistive magnetohy- 
drodynamic (MHD) models, which is parameterized solely by the dimensionless Lundquist 
number, 

^^fioLcsVA^ (1) 
V 

as a starting point of the discussion for magnetic reconnection. In Eq.Q, Lcs is the half 
length of the reconnecting current sheet, and can be taken as Lcs = where L is the 
plasma size and < e < 1/2 (the choices of e are discussed in Sec. VI. A.). Va is the Alfven 
velocity based on the reconnecting magnetic field component and t] is the plasma resistivity 
due to Coulomb collisions. The well-known Sweet-Parker modeP^ predicts reconnection 
rates as an explicit function of S*, 

^ - — (2) 

where Vr is the reconnection inflow speed. When collisions are sufficiently infrequent or 
5* is sufficiently large, physics beyond resistive MHD becomes cruciaP, leading to a fast 
reconnection rate nearly independent of 5*. A large body of the work in the past decades, 
therefore, has focused on reconnection either in collisional or collisionless limit as summarized 
by recent reviews'^'^. 

The collisional MHD description provides a good description of magnetic reconnection 
for plasmas in which all the resistive layers remain larger than the relevant ion kinetic scale. 
For example, without a guide field {i.e. anti-parallel reconnection), the transition between 
collisional and collisionless reconnection occur^^'^ when the current sheet half thickness 
predicted by the Sweet-Parker model approaches 



where di = c/upi is the ion skin depth. By properly varying both L (and thus Lcs) and 
7] (through changing e.g. electron temperature), S and di can be kept constant while the 
relative magnitude of 6sp to d^ can be reversed, leading to dramatic differences in the 
structure of the reconnection layer along with clear changes in the magnitude and scaling of 
the reconnection rate. This qualitative change can be characterized by the effective plasma 
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size which is defined by 

so that the second equahty in Eq.(|3| can be written as 

S = t^W (5) 

In the case of a finite guide field, the transition occur^n^H -vi^hen 5sp = Ps where ps = 
a/ (Tj + Tf^nii/ QiBx is ion sound gyroradius, Te and Tj are electron and ion temperatures, 
Bt is the total magnetic field including both the reconnecting and guide components, and 
and Qi are the ion mass and charge. In the case of anti-parallel reconnection with upstream 
plasma Pup <^ 1, ps will be equal to di by the virtue of the force balance across the current 
sheet, if the reconnecting magnetic field and the temperatures at the current sheet center 
are used to calculate ps- (We note that when P^p 3> 1, the transition scale for 5sp is less 
clear since Ps{— a/ f^updi) is separated from di.) Therefore, the boundary between collisional 
and collisionless reconnection is defined by Eq. ([s]) regardless of the presence of a guide field 
when the definition of plasma effective size is modified to 

A.^. (6) 

Ps 

Thus, the collisional and collisionless reconnection phases are distinguished in the parameter 
space of (A, S). This is illustrated as the black line in the phase diagram in Fig. [l] assuming 
e = 1/2. We note that the term "collisionless reconnection" is used in this paper for the 
reconnection process dominated by the effects beyond collisional MHD, such as two-fluid 
effects, ion and electron kinetic effects. Among these, the electron kinetic effects should 
become important in a similar parameter space deflned by the black line as shall be discussed 
in Sec.V. 



II. SINGLE AND MULTIPLE X-LINE COLLISIONAL RECONNECTION 

For plasmas larger than those specifled by Eq.([5]), it would appear that the collisional 
MHD description might be valid despite the large S. It shall become clear later in Sec. III. B., 
however, collisional models are not sufficient for describing reconnection in these regimes. 
Discussions on collisional reconnection have been long dominated by debates between the 
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FIG. 1: A phase diagram for magnetic reconnection in two dimensions. If either S or the nor- 
mahzed size, A, is small, reconnection with a single X-line occurs in collisional or in collisionless 
phases. When both S and A are sufficiently large, three new multiple X-line phases appear with 
magnetic islands. The dynamics of new current sheets between these islands are determined ei- 
ther by collisional physics or by collisionless physics (see Sec. Ill and Sec. IV.). The conditions for 
electron runaway are shown as red lines (see Sec. V). The locations for reconnection in Earth's 
magnetosphere, solar corona, solar chromosphere, and solar tachocline are also shown. The existing 
experiments, such as Magnetic Reconnection Experiment (MRX), do not have accesses to these 
new phases. A next generation reconnection experiment (NGRX) is required for such accesses to 
these new phases directly relevant to reconnection in heliophysical and astrophysical plasmas. 

Sweet-Parker modeP^ and the Petschek model-^, both of which, however, are unsatisfac- 
tory. The Sweet-Parker model has been verified numericallj^^^'^ and experimentallj^^ at 
relatively small values of S, but it predicts reconnection rates too slow to be consistent 
with observations of larger S plasmas. On the other hand, the Petschek model, invoking 
slow-mode shocks, predicts rates consistent with observations but it requires a localized re- 
sistivity enhancement in simulationJ^SHSS g^j^^ j^^^g been verified experimentally. The 
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origin of the localized resistivity enhancement is hypothesized to be kinetic in nature, but 
the underlying mechanisms still remain illusive. While signatures of slow-mode shocks have 
been reported in the Earth's distant magnetotaiP^, large-scale hybrid (kinetic ions and fluid 
electrons) simulations have revealed significant discrepancies in the expected structure of 
the discontinuities due to the strong ion temperature anisotropy that is naturally generated 
in these configur at ion^^^'^. 

There is a growing body of work that suggests it may be necessary to move beyond these 
steady-state models in order to understand the dynamics of magnetic reconnection in large- 
scale collisional plasmas. In particular, for sufficiently high Lundquist numbers, resistive 
MHD simulations feature highly elongated layers which breakup into muliple X-lines sepa- 
rated by magnetic islands (or plasmoids) ^^ * ^^* ^. These multiple-X line models are inherently 
time dependent, often generating impulsive reconnection consistent with observations such 
as Flux Transfer Events (FTE)!^. The plasmoid-like structures are also observed in Earth's 
distant magnetotail during substorm^ and in the current sheet during solar Coronal Mass 
Ejections (CME)^'^. Although the multiple-X line models were also applied to explain these 
observed plasmoids in the magnetotaiP^'^ or on the solar surfacd^^'^, they did not receive 
much attention until recent theorj^ and numerical simulationd^^^ offered detailed predic- 
tions concerning the break-up of Sweet-Parker layers to the so-called plasmoid instability, 
which produces numerous secondary magnetic islands. Although the time-averaged rates 
can be still differentpSESESl depending on the detailed divisions within (see Sec. Ill below), 
all of them are definitely much faster than the Sweet-Parker rate [Eq.(|2|]. Thus, the multiple 
X-line reconnection, associated with the plasmoid instability, constitutes a new reconnection 
phase within the collisional reconnection regime. Recent MHD simulations indicate that the 
critical Lundquist number, Sc, for the onset of the plasmoid instability is approximately 

Sc ~ lo^ (7) 

shown by the green line in Fig. [1} However, the precise value of Sc probably depends on 
the level of pre-existing fiuctuations in the plasmaP^^^'^^'^. We note that the green line is 
necessarily stopped at the low A end defined by Eq.(|5]) due to the invalidity of MHD models 
in the coUisionless phase. 
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FIG. 2: In the regimes of high Lundquist number and large plasma size, the plasmoid instability 
gives rise to a hierarchy of interacting current sheets and islands. The above sketch gives the 
notation used here for describing this hierarchy. 



III. MULTIPLE X-LINE COLLISIONAL AND COLLISIONLESS RECONNEC- 
TION 

Until quite recently, the boundary between collisional and collisionless reconnection was 
thought to be given by Eq.Q. This may not be true anymore when the current sheet is 
unstable to the plasmoid instability, forming thinner current sheets which may be further 
subject to new plasmoid instability leading to yet thinner current sheets in a hierarchical 
fashion as proposed by Shibata and TanumeP^. Eventually, these new current sheets can 
approach the ion kinetic scales triggering collisionless reconnection as recently demonstrated 
by full kinetic simulations with a Fokker-Planck treatment of Coulomb collision^. There- 
fore, the multiple X-line phase can be further divided into a phase involving only collisional 
physics (i.e. purely resistive MHD) and a phase involving both collisional and collisionless 
physics, which we denote as the multiple X-line hybrid phase. 

The boundary line in the (A, S) space between these two sub-phases depends on the de- 
tailed physics of unstable current sheets. The main uncertainty originates from the question 
of how many islands, on the average, remain within the unstable current sheet at any given 
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time. According to linear analytic theorj^^QES^ i^j^e number of secondary islands scales as 

N^{ff, (8) 

where a = 3/8. This linear prediction has been carefully verified in simulations designed to 
study the initial breakup of the Sweet-Parker layerd^^EU^ However, nonlinearly many more 
islands are observed in the simulationd^^ESEIM] corresponding to scaling parameters in the 
range a = 0.6 — i- 1. A possible explanation for this discrepancy is that, in these nonlinear 
simulations, the break-up of the original Sweet-Parker layer leads to new current sheets 
between the islands which are also unstable to the same plasmoid instability as illustrated 
in Fig. |2| and the islands on more than one level in the hierarchy were counted (see below 
for more discussions). 

At the present time, it appears that there are only two ways to terminate the downward 
progression in this hierarchy: (1) either the local Lundquist number of the new current 
sheets falls below the critical value for the plasmoid instability or (2) the new current sheets 
approach the ion ps scale where collisionless effects dominate. One can make quantitative 
predictions regarding these two possible outcomes with just a few simple assumptions. 



A. Multiple X-line collisional reconnection 

As illustrated in Fig. [2| we start by defining the half-length and thickness of the top 
level Sweet-Parker by Li = Lqs and Si = Ssp [Eq.([3])], corresponding to a macroscopic 
Lundquist number of 5*1 = 5 [Eq.Q]: Si = {Li/5iY. At this top level, the development 
of the plasmoid instability gives rise to A^^i = (Si/^c)" islands, which breaks the original 
layer into new sheets with length L2 = Li/Ni. We assume these new layers are governed by 
the Sweet-Parker scaling relationships and are also susceptible to the plasmoid instability 
in the same manner. Then, the number of islands generated within the second level of the 
hierarchy is A^2 = {S2/ Sc)" where S2 = (-^^2/^2)^ is the Lundquist number of the new sheets, 
assuming the same reconnection magnetic field strength upstream. Therefore, the jth level 
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quantities in the hierarchy are related to the (j-l)th level quantities by 

i = 4l = , ' = (11) 

where the number of islands in the jth level, Nj, is given by a recursion relation, 



Sc J \ J 

= iVji"") = ivji"")' = ... = 7V(^-")'"\ (12) 

If we terminate the hierarchy at the jth level, then the total number of islands in the system 
J\fj corresponds to the product of the islands in all the levels 

Afj = iViiVsiVs ...Nj 

^<^ = (|) (13) 



where 



ft.x:(i-«)"=i^ii^. (14) 

n=0 

Note that as the hierarchy becomes deeper j » 1 it converges towards ^oo — oT^ . We can 
now conveniently express the scaling for the total number of islands up through jth level in 
terms of the global Lundquist number 

while the Lundquist number of the new current sheets at the jth level is 



S 



- N,.iN,.,...Ni - N,^i ■ ^^^^ 

Notice that for a < 1 this result implies Sj > Sc for any finite level in the hierarchy, which 
implies the new levels are always unstable to the plasmoid instability! Strangely enough. 
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this implies the hierarchy has infinite depth for a < 1, and only terminates in the limit 
j — oo where Afoo = Si/Sc and S^o = Sc- However, our basic scaling assumption for the 
number of islands at each level Nj = {Sj/Sc)°' becomes invalid when Nj decreases towards 
unity. To make a reasonable estimate for the number of levels in the hierarchy, one should 
consider a cutoff Nj ~ Nmin below which continuous scaling arguments are meaningless. For 



example, by setting Nj > Nmin in Eq- ( 12 ) as a cutoff, the maximum level in the hierarchy 
is 

Jma. i+ In(l-a) ' ^^^^ 

which diverges logarithmically as A^mm 1- Picking some reasonable cutoff A^mm ~ 2 and 
assuming a = 3/8 will terminate the hierarchy fairly quickly {jmax = 2 — )■ 8) for most any 
conceivable Lundquist numbers (see Sec. VI). 



Therefore, the scaling of the total number of islands in the hierarchy [Eq. (15)] depends 



on the maximum number of levels, jmax [Eq. (17)], and the island number scaling power 
index from one level to the next level, a. Applying this estimate to the reported numerical 
MHD simulations^^ yields jmax ~ 3, using a = 3/8. This leads to the predicted scaling of 
■^jmax ~ ('S'i/5'c)''''^^. However, the linear scaling of a = 3/8 does not necessarily apply in the 
hierarchy model where the nonlinear evolution of islands at one level is required to generate 
new current sheets for the islands at the next level. Using a = 0.8, for example, leads to 
jmax ~ 2 and Mj^^^ ~ {Si/ S^)^-^^ . These scalings are not very far from the reported linear 
scaling of S*"^ given the large uncertainties that still exist (see Fig. 5 in Ref. 51). As the 
hiearchy becomes increasingly deep, the precise value of a no longer matters and the result 



approaches the linear scaling of S as evident from Eq.(15), consistent with earlier heuristic 
argument j ^^^ ' ^^ ' ^^ l 



B. Transition to multiple X-line collisionless reconnection 



There is a second way to terminate the downward progression in the hierarchy before 



reaching the maximum level estimated by Eq. (17). As discussed in Sec. I, this occurs when 



the thickness of a current sheet at a level j {> 1), given by Eq. (11), approaches the ion 
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sound radius to trigger coUisionless reconnection: 



Si 5i 



(18) 



Using Eq. (15), this can be expressed in terms of the critical global Lundquist number, 



S — as a function of plasma size, A, as 

^ = S'(^--)/2 (eA)^+% (19) 

where 



2 - (1 - 

which vanishes in the hmits of a = 1 or j = oo for < a < 1. The blue line in Fig. [T] 
shows the transition boundary in the (A, S) space from multiple X-line coUisional phase to 
multiple X-line hybrid phase in these limiting cases, 



S = V^eeA (20) 
where Sc = and e = 1/2, consistent with the previous heuristic argument d^^ ' ^^ l 
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For more realistic scenarios with a < 1 and a finite number of levels, the appxoimate 
transition boundary can also be estimated. For any given S{> Sc), the deepest level, jmax, in 



the hierarchy is given by Eq. (17). Then using this jmax, the maximum A for current sheets 



at the deepest level in the hiearchy to remain collisional is determined by Eq. ( 19 ). Examples 



for a = 3/8 and a = 0.8 are shown in Fig. [3] as a dotted line and a dashed line, respectively. 
The steps in these lines correspond to the increases in jmax, but overall they are only slightly 
above the limiting cases of c = 0. Lastly, we note that the physics across this boundary 
is vastly different: on the collisional side the reconnection is completely determined by 
collisional MHD physics while on the hybrid side, both collisional and coUsionless physics is 
important. This is consistent with the boundary determined by electron runaway conditions 
which also indicate the change in the required physics (see Sec.V later). The reconnection 

— 1/2 

rate, on the other hand, is given by the Sweet-Parker rate of ~ 0.01 on the collisional 

side while on the hybrid side, the reconnection is faster, but not by a large amount, at the 
collisionless rate of 0.01 — 0.1. 



IV. SINGLE AND MULTIPLE X-LINE COLLISIONLESS RECONNECTION 

The last part in our phase diagram concerns the fact that the single X-line current sheet 
in the collisionless phase [defined by Eq.(|5|] may be also subject to secondary collisionless 
tearing instability. Unlike the MHD counterpart, we are not aware of any analytic work in 
this area although there have been numerical demonstration^^SHSZl g^^^j some observational 
evidencd^. In the collisionless limit, sufficiently large kinetic simulations suggest!^ that the 
critical size for the secondary island formation in the extended current sheet as a result of 
nonlinear evolution is, 

A = Ac ~ 50, (21) 

which is shown as the orange line in Fig. [l] 

In principle, multiple X-line collisionless reconnection can also occur in the hybrid phase 
if the effective plasma size is larger than Ac at the hierarchy level when the current sheet 
thickness reaches ps- But it turns out that the condition for such transition is almost identical 



to Eq. (20) as shown below. The effective plasma size at the jth level in the hierarchy, Aj, 
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is given by 
Equating this to Ac yields 



A. = ^ = -T^ -4lV"'"°'" ■ (22) 



which reduces to S = {Sce/Xc)\ in the hmiting cases of a = 1 or j = oo. This condition 



is different from Eq. (20) by only a factor of 2 when Sc = 10^ and Ac = 50. Thus, the 
parameter space for the multiple X-line collisionless reconnection is very limited within the 
hybrid phase. 

V. ELECTRON RUNAWAY CONDITION 

In real plasmas, the whole concept of collisional resistivity (and thus Lundquist number) is 
restricted to parameter regimes in which the reconnection electric field is small in comparison 
to the Dreicer runaway limitP' given by 

En ^ (24) 
e 

where z/gj is electron-ion collision frequency. To put the phase diagram into better perspec- 
tives with regard to its previous versiorJ^, it is important to understand where runaway 
conditions are unavoidable and how these boundaries correlate with kinetic-scale transitions 
already discussed. For a Sweet-Parker reconnecting current sheet, the reconnection electric 
field is given by 

„ , TflpVpi Bfi 

ER = rij = ^^^^ 25 

where rj is classical resistivity, j is the peak current density, and Br is the reconnecting field 
component. Using the relation for the plasma /3 

\ 2 

ps 



/3 = 2(^) (26) 



where Tg = T, is assumed, we have 

Er 2^/2 Br irue Ps 



(27) 



Ed /3 Bt]I nii Ssp 
For the single X-line reconnection phase, Ssp is given by Eq.([3]), and Er/Ed = 1 leads to 

/32 f BtY x2 



128 \BrJ me ^ ' 



12 



where Bt is the total magnetic field. Assuming (3 = 0.01, Bt/Br = 10, rrii/me = 1836, the 
above equation becomes 

S = O.UX^ (29) 

which is below the boundary line defined by Eq.([5]) but within only a factor of 2 when 
e = 1/2. Therefore, the electron runaway condition is well coincident with the boundary 
line between coUisional and coUisionless phases. 



For the multiple X-line reconnection phase, we can use Eq.(18) with a = 1 or j = oo for 
simplicity, yielding 

5=^^,/^V^A. (30) 
Using the same example parameters as above, we have 

S = 0.38^/KX (31) 



which is again within a factor of 2 from Eq.(20) when c = and e = 1/2. Red lines in 



Fig. [T] illustrate the boundaries for electron runaway conditions, which separate collisional 
and coUisionless reconnection phases for both single X-line and multiple X-line geometries. 
The significance of the red lines in the phase diagram is that they separate the regime where 
reconnection can be described by collisional physics alone from the regime where coUisionless 
physics is required. The alignment of red lines with either the black line or blue line is 
consistent with the transitions from collisional reconnection to coUisionless reconnection, 
regardless whether it takes the form of the single X-line or multiple X-lines. 

Besides MHD models and fully kinetic models. Hall MHD models and hybrid models (fluid 
electrons and kinetic ions) are often used to study the transitions between collisional and 
coUisionless phase As demonstrated by the comparative studies between 
different model^. Hall MHD models and hybrid models can capture the qualitative bound- 
aries between these two phases. However, the coincidence of electron runaway conditions 
with the transition boundaries between collisional reconnection and coUisionless reconnec- 
tion raises questions on the suitability of these fluid models when they are used to study 
detailed dynamics near the transition j^ * ^^ * ^^ * ^ -^*^. The detailed electron kinetic dynamics 
become important in these regimes but are not yet accurately treated in fluid models. In 
particular, the transition into the multiple X-line hybrid phase unavoidably leads to runaway 
electric flelds {E > Ed) as illustrated by the red line in Fig. [l] Fully kinetic simulations 
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including the collision operator have demonstratecP^'^ that the mechanism breaking the 
frozen-in condition changes rapidly across this transition, from ordinary resistivity in the 
sub-Dreicer collisional limit {E <^ Ed) to off-diagonal terms in the electron pressure tensor 
for the runaway regime. Once this transition to runway electric fields has occurred, it is 
unlikely to be reversed as suggested by Hall MHD model^^*^ until fast reconnection has 
depleted the available flux. Indeed, large-scale collisional kinetic simulationd^^l have demon- 
strated that resulting electron layers in this runaway regime can become highly extended 
and are unstable to secondary magnetic islands in a manner similar to previous collisionless 
simulation^SSll] p^^^j-^er insights emerge when we divide 5* by A, yielding 

S ^ ^/^^ 

A 7] 

which is simply the Lundquist number based on p^. It has been suggested that there exists 
a critical value of S/A ~ 50 where the dynamics can revert from Hall dynamics (kinetic) 
back into the Sweet-Parker regim#'^. However, notice that Eq.(31) implies that electron 
runaway will occur for S/X > 40 beyond which simple resistivity models are known to 
break down. It remains an outstanding challenge to properly model this dynamics within 
two-fluid approaches, but comparative studies between these different models should be 
useful to provide guidance on reliable two-fluid models which can be practically used for the 
detailed investigations of the phase diagram at large S and A values. 



VI. DISCUSSION 



While the simple 5* — A diagram conveniently summarizes much of the present knowledge 
regarding the dynamics of magnetic reconnection, there are a variety of other factors that 
may significantly influence reconnection which have not yet been discussed. For example, 
we have largely avoided the onset question {i.e. how reconnection gets started) which is 
likely very different for the various parameter regimes. In particular, the detailed properties 
of tearing instabilities in various regions of parameter space may play some role in the onset 
phase of reconnection, but this subject is beyond the focus of the present paper. As another 
example, the structure of the large-scale initial condition can also influence the reconnection 
dynamics, including features such as asymmetric current layers and velocity shear. There 
has been some work on both of these issues, but many uncertainties remain. Below we 
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specifically discuss external drive dependence and three-dimensional effects, followed by the 
discussions on the locations of various plasmas in the phase diagram. 

A. Dependence of external drive 

One might consider that the reconnection process should depend on the external drive 
or on how much free energy available in the system for reconnection. This is especially 
relevant when reconnection is modeled by a local box around the diffusion region. In real 
applications, the boundary condition may significantly influence the reconnection process 
within. However, we point out that the definition of the Lundquist number given by Eq. ([T]) 
already takes into account of the available free energy in the system: the half length of the 
reconnecting current, Lcsi which is taken as a fraction of the system size, L, is used. If the 
system is completely relaxed without free energy for reconnection, then Lcs = even L can 
be very large. Having Lqs on the order of L implies that the available free energy is near its 
maximum. One can imagine a time evolution when the system is driven from its completely 
relaxed state with Lcs = to Lcs — L/2 to reach a maximum S", and the dependence on 
the drive is actually reflected in the magnitude of 5* already. If the free energy is less than 
its maximum in a given system, S should be less than its maximum value even L is still 
same. 

B. Influence of realistic three-dimensional dynamics 

The ideas leading to the phase diagram in Fig. [T] are largely based on two-dimensional 
(2D) models and simulations of reconnection. At present time, very little is known regarding 
how reconnection will proceed in large three-dimensional (3D) systems - either in fluid or 
kinetic parameter regimes. To begin with, the whole idea of magnetic islands relies upon a 
high degree of symmetry, which can be achieved in laboratory plasmas (and 2D simulations) 
but is unlikely to occur in space and astrophysical plasmas. Instead, extended flux ropes 
are the natural 3D extension of magnetic islands, and the manner in which these can form 
and interact is much more complicated than in 2D models. This is true of both the primary 
flux ropes which may form due to tearing instabilities, and also secondary flux ropes which 
can form in the new current sheets that arise from the nonlinear evolution of the primary 
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FIG. 4: Results from a kinetic simulatiorPSl of guide field reconnection showing the formation 
and interaction of flux ropes as illustrated by an isosurface of the particle density colored by 
the magnitude of the current density along with sample magnetic field lines (yellow) . Simulation 
parameters are rrii/me = 100 with the initial guide field equal to the reconnecting field. The 
domain size is TOdj x TOdj x 35(ij corresponding to 2048 x 2048 x 1024 cells and ~ 10^^ particles. 

flux ropes. The advent of petascale computers over the last few years is permitting 3D 
kinetic simulation^ to explore these ideas for guide field reconnection geometries, where 
tearing modes can be localized at resonant surfaces across the initial current sheet. These 
initial simulations together with Vlasov theory have demonstrated that the spectrum of 
oblique tearing modes within ion-scale layers is simpler than previously thought!^, but the 
resulting flux rope dynamics is still quite rich. Furthermore, the nonlinear development 
of the primary flux ropes produces intense electron-scale current sheets near the active x- 
lines and along the separatrices. As illustrated in Fig. |4| in 3D simulations these layers are 
unstable to the formation of secondary flux ropes over a broad range of oblique angles. The 
continual formation and interaction of these flux ropes gives rise to a turbulent evolution 
that is significantly different than 2D models. However, the full implication of these results 
will take years to sort out; researchers are just beginning to scratch the surface. 

In addition to these fundamental issues associated with island (flux rope) formation, there 
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are a wide range of other processes to consider in 3D which may potentially influence the dy- 
namics of magnetic reconnection. These processes include the lower-hybrid drift instability, 
driven by the strong diamagnetic currents, streaming instabilities, modes driven by either 
electron or ion velocity shear, and a range of kinetic instabilities driven by temperature 
anisotropy. Even in MHD, influence of a pre-exiting turbulence on reconnection remains an 
outstanding issue in both 2DESI and 3DEI studies. Huge challen ges remain in understanding 
the role these various process play in reconnection, and how they might change the phase 
diagram in Fig. [TJ One of the most long-standing ideas is that instabilities may modify 
the dissipation physics within electron-scale regions. However, there are other possibilities 
to consider including non-linear couplings between electron and ion-scale features, or the 
possibility that these instabilities may seed the formation of new flux ropes. 

C. Reconnection in heliophysical, astrophysical and laboratory plasmas 

Despite these rather serious caveats discussed in the previous section, it is interesting to 
place plasmas from laboratory, heliophysics and astrophysics in the phase diagram. In Fig. [T| 
some heliophysical and laboratory example plasmas are shown. In this section, results from 
a more extensive survey of astrophysical plasmas are summarized in Table [T] and Fig. |5] with 
references from which typical plasma parameters were taken. In general, these parameters 
are associated with large uncertainties due to limited measurements available from these dis- 
tant plasmas and crude models used for the estimation. Extreme astrophysical condition^, 
such as special relativity and radiation, are not taken into account here since these effects 
on coUisionless plasmoid instabilities is just beginning to be explored numerically *^"^ * ^"^ . On 
the log-scales as in Fig. [5} however, even an order of magnitude of the uncertainty does not 
change the location of these plasmas by much in the phase diagram. 

The cases shown in the phase diagram can be roughly grouped into three groups. The flrst 
group includes high temperature fusion plasmas and Earth's magnetospheric plasmas. These 
plasmas are completely in the coUisionless phases, either with a single X-line or multiple X- 
lines, depending on whether the plasma effective sizes. A, are larger than the critical Ac. The 
plasma for the Sgr A* flares may also belong to this group. 

The second group of plasmas cluster along the black line separating multiple X-line col- 
lisionless phases and multiple X-line hybrid phase. It spans over huge ranges from solar 
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TABLE I: Key parameters of various plasmas from laboratory, heliophysics and astrophysics. Un- 
less explicitly stated, assumptions are (1) e = 1/2, (2) the reconnecting field is 1/10 of total 
magnetic field, = O.IBt, (3) equal electron and ion temperatures, Tg = Ti, and (4) ions are 
protons. We note that there are opinions that the plasmas in Crab pulsar wind and radio lobes 
are nonthermal so that temperature may not be a good descriptiorf^^'^. There are some labora- 
tory experiments which are not listed: flux rope experiment^^^'^ with S* ~ A ~ 10^, and plasma 
merging experiment j'^^l^^l with S = 10^ — 10'^ and A = 10^ — 10^. 



location 


plasma 


size(m) 


r,(eV) 


?ie(m 


i3T(Telsa) 


5 


A 


Notes 






0.8 


10 


1 X 10"" 


0.1 


3 X 10'* 


1.5 X lO'^ 


f — 1 /4 T- — T /2 fto — "^Rt 






0.4 


25 


1.5 X 10"** 


0.044 


3 X 10^ 


4 X 10" 


e = l/A,Ti = 5eV, Ar+ 






2 X lO"'* 




5 X lO'^"'' 


100 


2 X 10" 


1 X 10" 




Lab 


mstEi] 


1 n 


1 3 X 10'* 


Q V 1 
y X lu 




O X lU 


6 2 X 10" 


±1 — OUUc V , , -L>/^ — xj.kjOjD'J^ 




TFTFp3 


0.9 


1.3 X 10" 


1 X 10^" 


5.6 


1 X 10** 


2.3 X 10^ 


Ti = 36keV, D+. Bn = O.OISt 




ITERpl 


4 


2 X lO" 


1 X lO'^" 


5.3 


6 X 10** 


5 X 10^ 


D+ Bo — n ni Bt 




NnT!)<|M] 


1.6 


15 


3 X 10"'' 


0.5 


2 X lO'"" 


1.2 X 10'* 


f — 1 /4 T- — T /2 Bd — n 2Bt 




N'l fig 11 G t O p 3,U.S 


V 1 

U X lU 


ouu 


1 X 10^ 


^ V 1 n^<^ 

O X lU 


U X lU 


Q V 1 n2 
y X lu 


IDR — JJt, VP.ZO/ ) 






6 X 10** 


600 


3 X lO'"" 


2 X IQ-* 


4 X 10^'* 


1.3 X 10'* 


Rr. — Hrj^ T ~ A 91:pV I'd '^'^'W 

J-J/^ — J-jJP J J- 1 — 'i.^JVc V , ^p.ZiOO ) 


Solar 


Solar WintPS 


2 X 10^" 


10 


7 X lO** 


7 X 10-" 


3 X 10^2 


2 X lO-' 


(p.92) 




Snlflr rnrnndS 


1 X 10"^ 


200 


1 X lO""' 


2 X IQ-'^ 


1 X IQi'* 


4 X 10^ 


IP- ' -^l 




Solfir CliromospliGr^^"^ 


1 X 10^ 


n ^ 
u.o 


1 X lO"^ 


2 X IQ-*^ 


1 X 10** 


V in** 

O X lU 


iiGutra,! particle effects are weali^^ ^ 




Solar TachoclinPiil 


1 X 10^ 


200 


1 X 102» 


1 


1 X 10" 


5 X 10"" 






Prntnstpllar Dislr'^iJ 


9 X 10^ 


3 X 10^^ 


6 X 10** 


2 X IQ-''' 


8 X 10'* 


1 X lO" 


J.J £^1 L\± L -L \J i ^ Cll ^*_J111S11_J11D 

iiiLiuuecp — ; -Lvig 




X-ray Binary DisksPEU 


4 X 10" 


75 


1 X 10^^ 


36 


3 X 10^ 


9 X 10** 


M = IOMq, L = 2h{R = lO^iJs), 






3 X 10" 


5 X 10"'' 


1 X lO'^" 


1 X lO" 


1 X 10"^ 


9 X 10^ 






CoronafPS 














{mplm.Q)T^. Tjcompton included^ 


Galaxy 


Crab Nebula FlarelSSMJ 


1 X 10" 


130 


10** 


10-7 


5 X 10^" 


2 X 10"" 


pair plasma, T from D'^/2jiQ = 2nr 






Gamma Ray BurstJ^U 


10" 


3 X lO-' 


2 X 10^5 


4 X 10" 


6 X 10"7 


2 X 10"f' 


pair plasma 




Magnetar FlarePEU 


10" 


5 X 10-' 


10"" 


2 X 10"i 


6 X 10"** 


5 X 10"' 


pair plasma, SGR 1806-20 




Sgr A* FlarefPEU 


2 X 10"" 


7 X 10** 


IQi-* 


10--* 


2 X 10^" 


5 X 10** 


L = 2R = 20Rs 




Molecular CloudfPlEI] 


3 X lO"'' 


10-3 


10" 


2 X 10-" 


1 X 10" 


7 X 10"2 


neutral particle effects includecP^ 
HCO+ 




Interstellar MedicPEU 


5 X 10"'* 


1 


10^ 


5 X 10-1" 


2 X 10^" 


1 X 10"" 


L=magnetic field scale height 




AGN DiskfPEZESJ 


2 X 10"" 


24 


8 X lO^-* 


0.5 


2 X IQi'* 


1 X 10"" 


M = IO^A/q, L = 2h{R = W^Rs), 
a = 10-^ M = lO'^^g/s 




AGN Disk Coronae^^ 


3 X 10"" 


5 X lO"' 


1 X 10"^ 


4 


1023 


3 X 10"" 


M = lO**Af0, R = Rs, Ti = 


Extra- 
















{mp/me)Te, i-jcompton includecp] 


galactic 


Radio LobeiPl 


3 X 10"9 


100 


1 


5 X 10-1" 


2 X lO^-' 


8 X 10"^ 






Extragalactic JetP^ 


3 X 10"9 


10" 


3 X 10"1( 


i 10-7 


6 X 10^" 


1 X 10"" 


3C 303 




Galaxy ClusterPS 


6 X 10"« 


5 X 10=* 


4 X 10" 


2 X 10-" 


2 X 10^-' 


6 X 10"" 


A1835 




log (k) 

FIG. 5: Various laboratory, heliophysical and astrophysical plasmas, in which magnetic reconnec- 
tion is believed to occur, are shown in the phase diagram. Other symbols are same as in Fig{T] See 
the text and Table H] for details. 



corona, accretion disk coronae, Crab nebula flare, to galaxy clusters, radio lobes, and extra- 
galactic jets. When 5* and A are both small, this same line separates single X-line coUisional 
phase and collisionless phase. With a single X line, the reconnection in the collisional phase 
was known to be much slower than in the collisionless phase. The plasma coUisionality was 
arguecP^ to regulate itself so that the plasma always stay near the marginal coUisionality, 
based on the reasoning that fast reconnection should effectively release magnetic energy 
evaporating nearby dense neutral gases (such as in the solar chromosphere) to increase den- 
sity and coUisionality until reconnection slows to a collisional rate. An alternative model was 
also proposed^ based on self-regulation of electron temperature to maintain marginal coUi- 
sionality through a similar but different reasoning: higher temperature lowers coUisionality 
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and fastens reconnection, and thus depletes quickly available magnetic energy and eventu- 
ally slows reconnection and cools the plasma while lower temperature increases collisionality 
and slows reconnection, and thus accumulate magnetic energy and eventually trigger faster 
reconnection and heat the plasma. 

However, at large 5* and A values for all plasmas in the second group, the marginality 
black line now separates multiple X-line collisionless phases and multiple X-line hybrid phase 
in the phase diagram. Now there is numerical evidenc^^^'^ that the reconnection rates in 
multiple X-line hybrid phase are as fast as the single X-line collisionless rate, consistent with 
the theoretical argument!^'' that the global reconnection rate is determined by a dominant 
reconnection site in the island hierarchy which should be collisionless in the hybrid phase. 
Therefore, the self regulation arguments for the collisionality mentioned above do not seem to 
hold at the large S and A phases, but much work still remains to be done. The accumulation 
of energetic particle populations is suggestecP^ as another player in the self regulation process 
of reconnection rate in the multiple X-line collisionless phase. Energetic particle populations 
should be regulated also by finite collisions in the hybrid phase, but detailed dynamical 
processes need to be investigated. 

The third and final group of plasmas shown in Fig. [5] occupy much of the multiple X- 
line coUisional phase: accretion disk interiors, solar chromosphere and tachocline, molecular 
clouds, gamma ray bursts and magnetar flares. It could be argued that they form a line 
slightly below but along the boundary (blue line) between the multiple X-line collisional and 
hybrid phases. It is conceivable that the self-regulation arguments for collisionalitjff^^'^ 
could be applied here since collisional reconnection dominates at the deepest level of the 
hierarchy on the one side of the boundary while collisionless reconnection dominates on the 
other side. In fact, it has been suggested through Hall MHD simulation^ that reconnection 
in the multiple X-line collisional phase is much slower than that in the hybrid phase although 
it is much faster than the single X-line collisional (Sweet-Parker) rate. However, the recon- 
nection rate is not so different: collisionless rates are around 0.1 and while the collisional 
rates are around I/a/S^ ~ 0.01 at the deepest level of the hierarchy. A key question here 
is what determines the overall reconnection rate in a hierarchy of islands and whether it is 
indeed dominated by the reconnection process at the deepest level^ or by the reconnection 
process at all levels in an integrated way. 

There are two special cases which do not belong to either of the above three groups: 
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protostellar disks and interstellar media. Protostellar disks have lowest 5* among the objects 
we surveyed and are located between the single X-line and multiple X-line coUisional phase. 
Interstellar media are right in the middle of the hybrid phase, and probably both collisional 
MHD physics and collisionless physics are important in charactering reconnection processes 
there ClS db part of the galactic dynamo. 

Lastly, we note that currently there are no laboratory experiments which can be used to 
study all of these new phases of magnetic reconnection. Laboratory experiments have been 
playing important roles in the reconnection research: confirming some leading theoretical 
or numerical models such as Sweet-Parkei'^ and collisionless reconnection model^^ while 
challenging others such as Petschek model; benchmarking state-of-the-art numerical simu- 
lationd^^nnSHQil (discovering 3D phenomensP^^'^^ studying flux rope dynamic j^^ * ^^ ', to name 
a few. As mentioned above, the main research tool on physics of new reconnection phases 
is numerical simulations using either full particle. Hall MHD, or resistivity MHD codes. 
Existing experiments, such as MRX, do not have accesses to these new phases which are 
important for the emerging themes of particle acceleration by magnetic reconnectiorP^'^. 
While numerical simulations, coupled closely with analytic theory, will continue to be a 
major player at this front, a next generation reconnection experiment (NGRX) based on 
the MRX concept is considered as a candidate for such a laboratory experiment!^. The 
parameter ranges for both MRX and NGRX are also indicated in the phase diagram for 
their coverages. 
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